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We try to calculate the interaction between vortices in the color-flavor locked (CFL) phase of 
Quantum Chromodynamics (QCD). The most fundamental strings in high density color super- 
conductivity are the non-abelian semi-super fluid strings. Using Abrikosov ansatz, one can show 
that the interaction between these vortices has a universal repulsion form when they are far from 
each other. The Ginzburg-Landau (GL) Lagrangian of CFL phase may be compared with the 
Lagrangian of Quantum Electrodynamics (QED) with three condensation states. The behavior of 
vortices of multi-component superconductivity, within the framework of GL theory was done for 
two condensation states (MgB2 superconductors) by variational method. The interaction between 
these vortices is attractive when they are far from each other and repulsive when they get close. 
We attempt to obtain the interaction between these vortices in such substances using numerical 
variational method for three condensation states (Fe-based superconductors). The ultimate aim of 
the research is to apply these calculations to the vortices in the color-flavor locked phase of QCD 
to be able to study the interaction between them. 
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Quark-quark pairing breaks the symmetry of QCD to CFL phase symmetry (G = SU {N)c x 
SU{N)f xU{1)q ^ H = SU{N)f+c X Zn)- Since the remaining symmetry is not abelain and 
the first fundamental group is not trivial, non-abelian vortices emerge [jl]] (G — )■ H ,H\{G /H) ^ 
/, TTi {G/H « U {N)) = Z). We are interested in the interaction between these non-abelian vortices. 
Abrikosov ansatz can be used when they are far enough. But the question is how these non-abelian 
vortices interact when they are close. By expanding the Lagrangian of the non-abelian vortices one 
can see that it is like the QED Ginzburg-Landau Lagrangian for multi-bands superconductors with 
a Josephson coupling. We plan to use the method of calculating the interaction between vortices 
introduced by Rebbi[Q]. Generalizing the G-L Lagrangian to the CFL phase of QCD, one gets 

^ = tr{D^y{D^) - mhr{^^'^) - Ai (trO^O)^ - Aitr [(<I>^<I>)^] - ^F^^-F"'^' , (1) 
where O is the matrix field corresponds to the pairing gap which its elements are defined by 

^ai^^aPr^ijkivJ'^CYsVl)- (2) 

These are the color and flavor antisymmetric representations of diquark Cooper pairs with spin 
zero, a, j3, and etc are flavor indices and /, j, and etc are color indices. C is the charge conjugation 
operator and Y indicates the quark field. The fundamental string, non-abelian string, is generated 
by both SU {N)c,f and the U{1)b generators. A string with a symmetry along the z-axis in polar 
coordinate (p, 0) can be shown by 



<I.(0,p) = exp (/f ) exp (^-/r;v2_i ^^^l^e) diag (/(p),g(p), • • • ,g(p)) = diag(^^^^^^ 



The interaction between vortices when they are close to each other has already been studied in 
QED|j^ |3|]. There exist two types of superconductors defined by the G-L parameter. It is type I 
if^ = A/<^<l/ \/2 . X is the London penetration depth and ^ is the correlation length. It is type II 
if k = X/^ > 1/ \/2. Possible interactions between QED vortices are attraction for superconductor 
type I and repulsion for superconductor type II. There is a possibility of attraction when the vortices 
are far from each other and a repulsion when they get close. This type of superconductor is called 
type 1.5, which is somehow new to the physics world. In these substances we have two condensa- 
tion states. Studying the behavior of substances with three condensations is interesting. We hope 
to use the results of this type of superconductors in the CFL phase of QCD. In this research, we use 
the variational method to study the interaction between vortices. First, we apply this method to the 
real superconductors. The Lagrangian of a substances with three condensations is 



iq^.l 2 I ft Im.l 4 I mi 



— Lr= 1,2,3 
+ (V X A)2 - 7i (»F*»F2 + 



iV-A)»F, 



(4) 



All coefficients and symbols are defined in 
couplings between condensations are equal (7 : 

na(i-) 



. For simplicity we consider the case when all 
7). We look for the solution = fi{r)e'"^ and 
A = ^^^Cq. The interaction will be realized if we know the solution of these equations for all 
range of distances. The following trial function is used 
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Figure 1: Left: Free energy density for a vortex with three condensation states. Right: Three condensation 
states and magnetic field versus distance from the center of the vortex. 

/i(r) = ^/T+ W + exp(- ^) ILo [fu r'/l\) , 

/;-2,3(r) = ^JH^+^ + ^W (-Tfe) iLo {fu r'/n) , (5) 

a(r) = l+exp(-^)i:;U («/ ^7^0' 
where f\j, f2j, /a,/ and a; are the variational parameters, rj is introduced in [||]. Using the 
Newton method with the iteration procedure vj''"^^'^ = v/"'^ [H^] ij^^j"'^' the variational pa- 
rameters are obtained. The superscript m represents the value at the mth step. D,- = j 5y. \y_y{m) , 
and the Hessian matrix is Hij = /dVidVj\ The stationary solution to the iteration pro- 

cedure corresponds to the minimum of the free energy. Figure |T] shows the free energy density 
versus the distance from the center of a voitex. The calculated free energy is finite, as expected. 
The plot on the right hand side of figure 1 shows the condensation states and the magnetic field of 
a vortex with distance. We choose <§i = 5 Inm, Ai = 25nm; ^2 = ^nm, A2 = 30nm and ^3 = 25nm, 
A3 = 5lnm so that the first condensation is of type I and the two others are of type II. At asymptotic 
distances from the center of the vortex, magnetic field vanishes and each condensation obtains its 
final value, depending on its correlation length and penetration depth. In the next step, we will put 
two vortices at different distances to calculate the interaction between them. Eventually, we would 
like to apply this method to the vortices of the CFL phase of QCD and study their interactions. 
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